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We study the sign problem of the fermion determinant at nonzero baryon chemical potential. For 
this purpose we apply a simple model derived from Quantum Chromodynamics, in the limit of 
large chemical potential and mass. For SU(2) color, there is no sign problem and the mean-field 
approximation is similar to data from the lattice. For SU(3) color the sign problem is unavoidable, 
even in a mean-field approximation. We apply a phase-reweighting method, combined with the 
mean-field approximation, to estimate thermodynamic quantities. We also investigate the mean- 
field free energy using a saddle-point approximation KM. 
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Model study of the sign problem in a mean-field approximation 



1. Introduction 

One of the current thrusts of hadronic physics is to understand the extreme conditions at high 
temperature and/or densities. The RHIC experiments reveal interesting features of a Quark-Gluon 
Plasma (QGP) phase above the critical temperature. In the core of a dense neutron star, various 
color superconductors should exist. Lattice calculations, based on Monte-Carlo simulation, are a 
powerful tool for the nonperturbative analysis of QCD. At zero baryon density, results from the 
lattice provide us with fundamental information, such as the phase transition temperature T c [0], 
the equation of state [Q], susceptibilities [Q], the behavior of correlation functions, and so on. 

On the other hand, at nonzero quark density lattice simulations have a serious sign problem: 
the quark determinant is complex in the presence of the baryon chemical potential, so Boltzmann 
weights are complex, and importantance sampling fails. There are ways to overcome the problem, 
including reweighting [||], Taylor expansion in the chemical potential and analytical continua- 
tion from imaginary values of the chemical potential [|7|]. These methods are applicable when the 
chemical potential is small, and the temperature high. The sign problem is even more intractable 
at low temperature and high density. In this work we study the sign problem in a toy model using 
a mean-field approximation. 

Let us first see how the sign problem arises. The quark determinant is 

det^(n q )=det[Y^ + Y 4 n q + m q ], (1.1) 

where D* 1 = d^ 1 — igA^ is the covariant derivative, m q is the quark mass and \L q is the quark chem- 
ical potential. The quark determinant is complex except for \L q = 0, due to a lack of gamma-five 
Hermiticity of Jt\\i q \. det^#(ju 9 ) = dety5^(/i 9 )y5 = {&&tj%(—ji q )}*. In and of itself, a com- 
plex quark determinant is not necessarily fatal. While the quark determinant is complex for any 
given Afx, the functional integral over is real for real observables. This is seen from the relation 

det[y M (<^ - ig{A») c + y^ q + m q )} = {det[y^ + y 4 ju, +m q }}*. (1.2) 

The real part of the quark determinant is C-even, and the imaginary part, C-odd. For a C-even (C- 
odd) observable, then, the imaginary (real) part of the determinant vanishes after integration over 
A^ . Accordingly the real problem is that the contribution of the quark determinant changes sign, 
depending upon A^ , and there is no known method to replace importance sampling. 

2. Model 

We analyze a simple model to see the sign problem in the mean-field approximation. The 
model is obtained by taking double limit of heavy mass, m — > °°, and large chemical potential, 
/i — > oo, keeping the ratio e = (e^' a /2m q a) Nt fixed. In the heavy quark limit all excited quarks 
are static, while antiquarks are suppressed at nonzero quark density, so that in the end, the quark 
determinant can be rewritten in terms of the Polyakov loop, 

e- S([L] = det [y^ + y 4 ^ + m 9 ] -> [det(l + eL)] N ^\ (2.1) 
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where L{x) = Ylx 4 L^ 4 (x) is the Wilson line, which is a color matrix. Here we set Nf = 4 in order to 
avoid the rooting problem with staggered quarks. The determinant can be explicitly calculated: 

Y\{\ + e 2 + 2ei) forSU(2) 
det^= { jj (1 + e 3 + 3 ^ + 3e 2 r) forSU(3) > W 

where the Polyakov loop is the trace of the Wilson line, £(x ) = trL(x) /N c in the fundamental repre- 
sentation. One can easily check that SU(2) color does not have the sign problem since the Polyakov 
loop is always real, — 1 < I < 1. For SU(3) case the Polyakov loop I is complex valued, and the 
determinant is complex. We analyze these two cases in the next section. 

For gluons, we take a simple action with nearest neighbor interactions between Polyakov 
loops, 

S g = -N^£(x)t(y). (2.3) 

n.n. 

/ is a parameter which can be interpreted as the temperature of the system. In the strong coupling 
expansion, J is related to the true temperture, T, through / = exp[— aa/T], where a is the string 
tension. It is known that this action reproduces the gross features of the phase transition without 
quarks; i.e., a second-order phase transition for SU(2) color, and a first-order phase transition for 
SU(3) color. In this work we leave J as a free parameter. 

3. Mean-field approximation 

At nonzero temperature the free energy is related to the functional integral as 

e-W =Z(P,n)=J DLexp(-S), (3.1) 

j3 = 1 /T. Assuming that the action S is real, exp(— S) is positive semidefinite, and in a mean-field 
approximation the free energy is: 



e -pvf 



= J DLexp(-S m f(x)+S m f(x)-S) = Z mf (j8)(exp(5 mf (x) -S)) mf 

> Z mf (j3)exp((5 mf (x) -S)mt) = e-WM. (3.2) 



The average (• • • ) m f is taken with respect to mean-field action S m f(x); x is a parameter of mean- 
field theory. In eq. (|0|), the first line is an identity, while the second line follows from Jensen's 
inequality, ((exp ff)) > exp((<^)). The mean-field free energy is larger or equal than the exact free 
energy for any x. The inequality ensures that the f(x) has a minimum at x = xq, 



df(x) 



= 0, and 

X=XQ 



d 2 f(x) 



> . (3.3) 

x=xo 



Our ansatz for the mean-field action is 



S m AL]^- X -l}l(x) + f(x)}. (3.4) 
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Once 5 m f [L] is known, and x determined, the expectation value of any observable 0[L] is given by 
integrating with repsect to the Wilson line, L, over the group measure, with the action of mean- 
field theory, S m f[L]: (0[L\) — {&[L]) m f. If the action is not real, inequality is not ensured, and 
convexity is violated. This how the sign problem manifests itself in a mean-field approximation, 
and occurs for three or more colors. Charge-conjugation symmetry is violated at non-zero quark 
density, which is seen from {£) ^ (f ). Here we note that both {£) and (£*) are real valued, from the 
argument in Sec.|l|. This difference has been observed in both lattice simulations |Q] and in other 
models [JSj] . It is necessary to extend the mean-field ansatz (3.4) to include two variables, x and y, 
in order to represent the difference between (£) and (£*): 

Smf = -\ %W) + £*(x)] -lZ[£(x) - f(x)\. (3.5) 

X X 

While the mean-field action is complex, x and y are real, so that after integrating over L, the free 
energy f m f(x,y) is a real function of x and y. Their values are then determined by requiring that the 
free energy is a stationary point. At \i q ^ 0, y ^ 0. It turns out that about the stationary point, while 
the free energy f m f(x,y) is minimal in the x direction, it is maximal in the y direction. That is, the 
solution is a saddle-point in x and y, consistent with Ref. [H]. 

The phase reweighting method is one way to deal with the complexity of the action. The 
magnitude of the quark determinant is C-even, while its phase is C-odd. Accordingly, the quark 
action is 

St = Sf ag + iQ[L], (3.6) 

where 

S mag = -£ln|l + £ 3 + 3££ + 3£f |, (3.7) 
&[L) = -£arg(l + £ 3 + 3££ + 3£ 2 f ). (3.8) 

With these definitions the expectation value of ff[L] is 

(0[L]) ~ (0[L}e- @[L] U/(e- i@[L] U. (3.9) 

Here 5 m f [L] or x is fixed from the free energy with the action S g [L] + 5™ ag [L] , so that x encompasses 
the information of S™ s [L] implicitly. This scheme is the same as what has been adopted in the lat- 
tice simulations of Ref. [B]. We compare these two methods for SU(3) color, and find qualitatively 
similar behavior. 

4. Results 
4.1 SU(2) 

We first consider SU(2) color, to see how a mean-field approximation works when there is 
no sign problem. We look for a phase transition by considering how the Polyakov loop changes 
as / increases. In the pure glue theory, the deconfining phase transition is known to be of second 
order for two colors, in the universality class of the Ising model. In our model, at £ = there is a 
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Figure 1: For SU(2), comparison of the model to lattice data, Fig. 2 of Ref. Left: the Polyakov 
loop versus the temperature parameter. Center: the Polyakov loop versus the density parameter. Right: the 
number density versus the density parameter. 



continuous transition at J = J c ~ 0.083, as indicated by the solid curve in the left figure of Fig. ||. 
The presence of dynamical quarks acts on the Polyakov loop variable as an external field which 
breaks the center symmetry. In fact, the results at e ^ in the left figure of Fig. [I] indicate not a 
true phase transition, but only crossover. 

Our mean-field outputs are to be compared with the lattice simulations in Ref. ^ : the center 
and left figures in Fig. [T] correspond to Figs. 1 and 2 of Ref. [^], respectively. We cannot expect 
exact agreement, because our ansatz for the pure gluonic action S g [L] is only a crude approximation 
of QCD, and in any case, we neglect the renormalization of the Polyakov loop in a mean-field 
analysis. Nevertheless, the agreement turns out to be surprisingly good, beyond naive expectation, 
if the parameter J is treated as an adjustable parameter as a fitting parameter. In this way, we fix 
J = 0.0042 and J = 0.04 to reproduce the SU(2) Polyakov loop only at e = 1 for 4/g 2 = 2.0 and 
4/g 2 = 1.5, respectively. We stress that we do not use the data of the Polyakov loop at e ^ 1, nor 
the results on the number density. Nevertheless, as clearly seen from Fig. [jj our numerical results 
fit all of the lattice data remarkably well. We conclude from this that the main corrections to our 
ansatz (2.3) can be represented by a shift in the parameter /. This gives us confidence in using a 
mean-field approximation for this problem. 

4.2 SU(3) 

We next consider SU(3) color. The Polyakov loop is compared in the phase reweighting 
method, and the saddle point approximation, in the left figure in Fig. ||; this is to be compared 
with the lattice results from Fig. 7 of Ref. [§]. We find a first-order phase transition for e = 
at J = J c = 0.132 and for e = 0.1 at J = J c = 0.123. Nonzero e smears the transition, so that it 
eventually ceases to be of first-order. The line of first-order transitions ends with a second order 
transition, which is then a critical end-point. For larger values of e there is only crossover. The 
global picture is consistent with results from a Potts model. In the left figure in Fig. |2| one sees that 
both reweighting, and the saddle-point approximation, have qualitatively the same behavior for the 
expectation value of the Polyakov loop. 

At ii q ^ 0, (£) ^ (£*). The observable (£ — £*) is C-odd, where the imaginary part of the 
fermion determinant is responsible for this difference. In the center figure of Fig. |2| we present our 
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Figure 2: For SU(3), comparison of the model to lattice data, Fig. 2 of Ref. [|]]. Left: comparison 
of the Polyakov loop between phase reweighting and the saddle point approximation. Center: difference 
between (£) and (£*), versus the temperature parameter. Right: expectation value of the phase of the quark 
determinant versus the density parameter. 



numerical results for the difference (£) — (£*} as a function of /. The difference is trivially zero 
at £ = and e = 1 where the fermion determinant is real. As long as the density parameter stays 
smaller than e ~ 0.5, a larger density parameter e leads to a bigger difference. For example, at 
e = 0.5 we find (£} - (£*} = -0.076, which is comparable to (£) = 0.073. 

One can intuitively understand why (£*) is greater than (£) at nonzero \l q , as seen in [Qj. It is 
because at nonzero quark density, the presence of quarks enhances the screening of antiquarks, so 
that an antiquark costs less energy [§]. 

Finally we present the results for the expectation value of the phase factor of the quark deter- 
minant, e~' @ . We plot (e~ l9 ) as a function of e in Fig. ^, where is the phase at each lattice site; 
d = -arg(l + £ 3 + 3£^ + 3£ 2 r), i.e. = 1^0. Comparing it with Fig. 9 in Ref . [g], we see that 
our results qualitatively reproduce the lattice data. For more quantitative agreement, we approxi- 
mate the phase factor by (e~' e ) 216 , taking the lattice volume of 6 3 = 216 from [9]. For instance, our 
7 = result has a minimum at e = 0.61 where {e~ w ) ~ 0.977, while we obtain 0.977 216 = 0.0066. 
The minimum value in Fig. 9 of Ref. [^] is 0.01, which is close to our value. 



5. Summary 

We have explored a simple model applicable in the limit of heavy quark mass and large chem- 
ical potential, and seen how the sign problem, at nonzero quark density, manifests itself in a mean- 
field approximation. All results from mean-field theory are reasonable, and are in quantitative 
agreement with lattice data. There is no sign problem for SU(2), and we find that a mean-field 
approximation works well, for both the quark number density and the Polyakov loop. 

For SU(3), we compared two methods, a saddle point approximation with a complex action, 
and phase reweighting. We find that both methods give qualitatively the same behavior for the 
expectation value of the Polyakov loop. The complex action implies that (£) ^ (£*}, which is 
also seen with phase reweighting. We computed (£) — {£*), as a function of the parameters for 
temperature and density. 
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While it may appear odd to find a saddle point in the mean field approximation, it is known 
that this happens in other theories, such as for a nonlinear sigma model. There, the constraint of 
the nonlinear model is eliminated by introducing a new field. The effective action, including the 
constraint field, is complex, so that the stationary points thereof are true saddle points. For QCD, 
the Aq field is a constraint field which imposes Gauss' law. This may provide a clue to resolving 
the sign problem with dense quarks. 
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